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ABSTRACT 

Estimates of the metagalactic hydrogen ionization rate from the Lya forest opacity 
in QSO absorption spectra depend on the complex density distribution of neutral hy- 
drogen along the line-of-sight. We use a large suite of high resolution hydrodynamical 
simulations to investigate in detail the dependence of such estimates on physical and 
numerical parameters in the context of ACDM models. Adopting fiducial values for 
cosmological parameters together with published values of the temperature of the IGM 
and the effective optical depth, the metagalactic ionization rates which reproduce the 
Lya effective optical depth at z = [2, 3, 4] are r H i = [1.3±g;|, 0.9±0.3, 1.0±g j] x 1CT 12 
s _1 , respectively. The errors include estimates of uncertainties in the relevant physi- 
cal parameters and the numerical accuracy of the simulations. We find the errors are 
dominated by the uncertainty in the temperature of the low-density IGM. The esti- 
mated metagalactic hydrogen ionization rate for the neutral hydrogen distribution in 
the current concordance ACDM model is more than four times the value inferred for 
that in an Einstein-de Sitter model of the same r.m.s. density fluctuation amplitude 
as. The estimated ionization rate is also more than double that expected from updated 
estimates of the emissivity of observed QSOs alone. A substantial contribution from 
galaxies appears to be required at all redshifts. 

Key words: methods: numerical - hydrodynamics - intergalactic medium - diffuse 
radiation - quasars: absorption lines. 



1 INTRODUCTION 

In recent years, high resolution, high signal-to-noise quasi- 
stellar object (QSO) spectra (e.g. Lu et al. 1996; Kirkman 
& Tytler 1997; Kim et al. 2002), combined with hydrody- 
namical simulations <CeJl992tlZhang. Anninos fc NormarJ 



119951: iHernauist et alJ Il99tj ; ilVIiralda Escude et alJ Il996t 
iTheuns et alll99Sft . have provided fresh insight into the na- 
ture of the intergalactic medium (IGM). It is now widely 
accepted that the Lya forest arises from absorption by den- 
sity fluctuations imposed on the warm photoionized IGM 
as suggested by analytical calculations (Bi, Borner & Chu 
1992, Bi & Davidson 1997), and is a natural consequence 
of hierarchical structure formation within CDM-like models 
(see Rauch 1998 for a review). 

One notable inconsistency with this paradigm has been 
the rather low values of Thi inferred from early hydrody- 
namical simulations of the Lya forest if scaled to appro- 
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priate assumptions for the IGM temperature (e.g. Rauch 
et al. 1997). However, the matter density S7 m , baryon den- 
sity Sib, amplitude of mass fluctuations as and the parame- 
terised Hubble constant h — Ho/WO kms _1 Mpc~ 1 , as well 
as the gas temperature and the Lya effective optical depth, 
t cff = — ln(_F) bs, all influence such an estimate. This has 
led to a wide range of estimates for Thi inferred from the 
Lya forest opacity using numerical simulations with differ- 
ent assumptions (Rauch et al. 1997; Theuns et al. 1998; 
Haehnelt et al. 2001; McDonald & Miralda Escude 2001; Hui 
et al. 2002; Schaye et al. 2003; Sokasian, Abel & Hernquist 
2003; Meiksin & White 2003, 2004; Tytler et al. 2004). This 
complicates the comparison with determinations of the ion- 
ization rate from estimates of the integrated emission from 
observed QSOs and/or galaxies (e.g. Haardt & Madau 1996, 
hereafter HM96; Madau, Haardt & Rees 1999) and estimates 
using the proximity effect (e.g. Bajtlik, Duncan & Ostriker 
1988, Bechtold 1994, Cooke, Espey & Carswell 1997, Scott 
et al. 2000) whic h both also have ra ther large uncertain- 
ties. Furthermore. ISteidel et alJ 11200 If) , in a study of Lyman 
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break galaxies at z ~ 3.4, suggested that the intensity of the 
ionizing background may be a factor of a few larger than in 
the model of HM96 due to a large contribution from star- 
forming galaxies. 

In this paper we present new estimates of the meta- 
galactic hydrogen ionization rate using Lya forest opacity 
measurements taken from the literature. Our approach ex- 
pands on previous studies in that we examine the depen- 
dence of Fhi upon the various cosmological and astrophys- 
ical parameters using a large suite of hydrodynamical sim- 
ulations. In section [5] we present an overview of our numer- 
ical code and our method for obtaining Fhi from artificial 
Lya spectra. We examine how Fhi depends on a wide vari- 
ety of parameters in sectional and compare our results to 
previous estimates using the same technique in section 0] In 
section|K|we present our best values for T_i2 at z = [2, 3, 4] 
and consider the implications for sources of ionizing radia- 
tion. Section E] contains our summary and conclusions. 



2 FROM MEAN FLUX TO IONIZATION RATE 

2.1 The fluctuating Gunn-Peterson approximation 

At redshifts where the Lya forest is commonly ob- 
served, density fluctuations in the low density, unshocked 
baryons closely foll ow those of the underlying dark matter 
llZhang et al]ll99sl) . The balance between photoheating of 
unshocked gas and cooling due to the adiabatic expansion 
of the Universe results in a power law temperature-density 
relation for the nor malised gas density A = p/{p) < 10, 
where T = TbA 7-1 llHui fc GnedirJll997ft . This is often re- 
ferred to as the effective equation of state of the IGM, for 
which typically, 10 3 ' 7 K < T < 10 43 K and 1.0 < 7 < 1.6. 
Additionally, in photoionization equilibrium the Lya optical 
depth, r, is proportional to the density of neutral hydrogen, 
rim- Neglecting collisional ionization, which is an acceptable 
approximation at the low densities associated with Lya for- 
est absorbers, nm oc /9 2 T _0 ' 7 /Fhi, where p is the gas density, 
T is the temperature and Fhi is the hydrogen ionization rate 
per atom. Combining these relations gives a simple power 
law expression relating the optical depth of Lya absorbers to 
the p hysical properties of the hydrogen gas llWeinberg et alJ 
1999). This relation has been termed the Fluctuating Gunn- 
Peterson Approximation (FGPA). 

The optical depth due to Lya absorption for a ho- 
mogeneous distribution of hydrogen g as in an isotropi- 
cally e xpanding Universe was discussed bv lGunn fc Peterson! 
(1965). The FGPA applies to an inhomogeneous IGM in 
photoionization equilibrium. Ignoring the effect of pecu- 
liar velocities and thermal broadening, the optical depth to 
Lya scattering is given by (e.g. Weinberg et al. 1999, Mc- 
Donald & Miralda-Escude 2001): 

(i + z) 6 (n b h 2 f p 
r = ro ro.^)r_ 12 (,) A (1) 

where is the baryonic matter density as a fraction of the 
critical density, h — iio/100 kms _1 Mpc _1 for the present 
day Hubble constant Ho, T is the gas temperature at nor- 
malised density A, r_i2(z) = Fhi(z)/10 -12 s _1 at redshift 
2 and ro is a constant. The index /3 = 2 — 0.7(7 — 1) is 
determined by the slope of the effective equation of state, 
and asymptotically approaches ~ 1.6 after reionization. For 



a flat Universe at redshifts z > 2 the Hubble parameter can 
be approximated by H(z) ~ _f/ofim 2 (l + z) 3//2 . Thus, given 
the correct density distribution, velocity distribution and ef- 
fective equation of state at redshift z, the Lya optical depth 
should scale according to the parameter combination: 



2 i,3T,-0.7 n -0.5 r -l 

p = il h h T S2 m T_ 12 



(2) 

for redshifts 2 > 2 (Rauch et al. 1997). This assumes that 
thermal broadening and peculiar velocities will not have a 
large effect on the optical depth distribution. At the den- 
sities and physical scales of the Lya forest absorbers Hub- 
ble b roadening is expected to be dominant llWeinberg et alJ 
Il997l) which should make this a reasonable assumption. 

We match the mean flux of artificial Lya spectra con- 
structed from hydrodynamical simulations to observed val- 
ues by linearly rescaling the simulated optical depths. If the 
FGPA approximation holds this corresponds to changing the 
parameter combination p by the same factor. We have deter- 
mined p in this way. If independent estimates are available 
for Qh, f2m, h, T and r e ff , the magnitude of T_i2 required 
to reproduce the observed mean flux can be derived. This 
assumes that the effect of absorption by fully neutral gas 
on the mean flux level is small, which should again be a 
reasonable assumption. 



2.2 Numerical code 

We proceed to investigate how F_i2 depends on several 
cosmological and astrophysical parameters using a suite of 
hydrodynamical simulations. The simulations were run us- 
ing a new version of the parallel TreeSPH code GADGET 
JSpringel. Yoshida fc Whitell200ll) . GADGET-2 has the op- 
tion of using a TreePM mode to perform the long range 
gravitational computations, which results in a significant re- 
duction in computational run-time. 

The conservat i ve entropy formulation of 
ISpringel fc Hernauistl <l2002l) is used to model the SPH 
particles, and radiative cooling is include d usin g the 
prescription of iKatz. Weinberg fc Hernauistl 1^996;). We 
have used the UV background specified by HM96 based 
on observed counts of QSOs, but with artificially increased 
heating rates to explore the effect of gas temperature on 
r_i2. Our fiducial model uses a heating rate a factor of 

3.3 greater than the HM96 value in the optically thin limit 
which resul ts in temper a tures close to those observed in 
the IGM dSchave et al.l l2000l : iBrvan fc Machacekl l2000l : 
iRicotti et alJ l200cl7 ! The fiducial temperatures at mean 
gas density we shall assume are T = [11200, 17800, 12500] 
K at z — [2,3,4], based on the results of ISchave et alJ 
(2000) which are consistent with He 11 reionization at 
2 ~ 3. At present there are only weak observational 
constraints on the index 7 of the t emperature density 
relation and i ts dependence on redsh i ft l|Ricotti et al. 2000; 
ISchave et all 120001 iMcDonald et alJ 1200 J) . The plausible 
range is 1 < 7 < 1.6 and in the following we will adopt 
7 = 1.3 as the fiducial value. Unless otherwise noted the 
temperatures of the gas have been rescaled such that the 
temperature at mean gas density To has the above fiducial 
values and the index of the power-law temperature-density 
relation for low gas densities is 7 = 1.3. At higher gas 
densities the temperature-density relation asymptotes to 
T ~ 10 4 K as atomic cooling becomes important. 
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Star formation is included using a simplified pre- 
scription to maximise the speed of the simulations. 
All gas particles with A > 10 3 and T < 10 5 K 
are converted into collisionless stars. Comparisons with 
identical simulations using the multi-phase model of 
ISpringel fc Hern auist] l|2003|) s how little differe nce in the 
Lva absorption lyieLIjaehneftTfe Springell20 04'l . Addition- 
ally, all GADGET-2 feedback options have been disabled. 
The effect of feedback from galactic winds on the flux distri- 
bution is uncertain but is believed to be small (e.g. Theuns 
et al. 2002; Bruscoli et al. 2003; Desjacques et al. 2004). Sim- 
ulations run by one of us (VS) indicate that the inclusion of 
winds only marginally alters the opacity of the Lya forest, 
such that the mean transmitted Lya flux predicted by dif- 
ferent models varies to less than half a per cent. There also 
appears to be no monotonic trend of the mean transmit- 
ted flux with wind strength. It is, however, conceivable that 
more realistic implementations of galactic winds within our 
simulations may lead to a multi-phase medium with dense 
cool clouds, which may increase the Lya forest opacity and 
thus the estimated hydrogen ionization rate. We neglect this 
possibility in this work. Note that Theuns et al. (2002) find 
the inclusion of winds alters the mean transmission by 10 
per cent. 

The simulation parameters are listed in Table □ For 
each model we have run a simulation with 200 3 dark matter 
particles and 200 3 gas particles within a lbh^ 1 comoving 
Mpc box with periodic boundary conditions. We find this to 
be the best compromise between accuracy and speed for our 
parameter study. For our fiducial model we have also run 
simulations with different box sizes and mass resolutions to 
test numerical convergence (see section l2~4l . These are listed 
in Table Note that the simulations with 200 3 particles 
used for the parameter study do not give fully converged 
estimates of T_i2 and we will later show results corrected for 
box size and resolution. In addition, we have run ACDM and 
SCDM simulations with 2 x 64 3 particles and box sizes of 
lO/i and ll.l/i comoving Mpc respectively, to mimic the 
simulations used in the seminal study of Rauch et al. (1997). 
The parameters are also listed in Tabled 

The simulations were started at z — 99, with initial 
conditions gener a ted us ing the CDM transfer functions of 
lEisenstein fc Hul lll999l) . The initial gas temperature was 
T = 227 K and the gravitational softening length was 2.5/i -1 
comoving kpc. The simulations were run on COSMOS, a 152 
Gb shared memory Altix 3700 with 152 1.3 GHz Itanium2 
processors hosted at the Department of Applied Mathemat- 
ics and Theoretical Physics (Cambridge). 

2.3 Method for determining T 12 from simulations 

Artificial Lya spectra are constructed at z — [2, 3, 4] from 
1024 random lines-of-sight (LOS) parallel to the box bound- 
aries (e.g. Theuns et al. 1998). Each LOS consists of 1024 
pixels. We mimic the noise on high resolution QSO spectra 
by adding Gaussian distributed noise with a total S/N = 30. 
The simulated optical depths we use are divided by a con- 
stant factor of 1.225 to correct the hydrogen recombination 
coefficient am used in our simulations, which is taken from 
ICenl l)l992t) . As noted by Rauch et al., t his value of am is 
around 20% larger than the calculation of Abel et al. ( 1997) 
within the relevant temperature range for unsaturated ab- 



Table 1 . Parameters used in our suite of hydrodynamical simu- 
lations where we assume a flat Universe with Q\ = 1 — f2 m . The 
last column lists the photoheating rate as a factor of the HM96 
value. 



Name 




n h h 2 


h 


CT8 


n 


uv 


Fiducial 


0.26 


0.0240 


0.72 


0.85 


0.95 


3.3 


Tl 


0.26 


0.0240 


0.72 


0.85 


0.95 


0.3 


T2 


0.26 


0.0240 


0.72 


0.85 


0.95 


1.0 


T3 


0.26 


0.0240 


0.72 


0.85 


0.95 


2.0 


Tl 


0.26 


0.0240 


0.72 


0.85 


0.95 


4.0 


Ml 


0.17 


0.0240 


0.72 


0.85 


0.95 


3.3 


M2 


0.40 


0.0240 


0.72 


0.85 


0.95 


3.3 


M3 


0.70 


0.0240 


0.72 


0.85 


0.95 


3.3 


M4 


1.00 


0.0240 


0.72 


0.85 


0.95 


3.3 


SI 


0.26 


0.0240 


0.72 


0.50 


0.95 


3.3 


S2 


0.26 


0.0240 


0.72 


0.70 


0.95 


3.3 


S3 


0.26 


0.0240 


0.72 


1.00 


0.95 


3.3 


S4 


0.26 


0.0240 


0.72 


1.20 


0.95 


3.3 


ACDM 
SCDM 


0.40 
1.00 


0.0150 
0.0125 


0.65 
0.50 


0.79 
0.70 


0.95 
1.00 


3.3 
1.0 



Table 2. Resolution and box size of our six additional simulations 
which have the same parameters as the 15-200 model used for 
our parameter study. The mass resolution used for our parameter 
study is also listed for comparison. 



Name 


Box size 


SPH 


Mass resolution 




comoving Mpc 


particles 


ft M0/gas particle 


15-400 


15k- 1 


400 3 


6.78 x 10 5 


30-400 


30ft" 1 


400 3 


5.42 x 10 6 


15-100 


15k- 1 


100 3 


4.34 x 10 7 


30-200 


30ft" 1 


200 3 


4.34 x 10 7 


60-400 


60ft" 1 


400 3 


4.34 x 10 7 


30-100 


30ft- 1 


100 3 


3.47 x 10 8 


15-200 


15ft" 1 


200 3 


5.42 x 10 6 



sorption in the Lya forest. Abel et al. base their fit for am 
on the data of Ferland et al. (1992). 

To determine F_i2 we rescale the optical depths of the 
artificial spectra to match the mean observed flux of the 
Lya forest at 2 = [2,3,4]. W e use the central v alues of 
r G fr from the fitting formula of ISchave et al] (l2003fl . based 
on a sample of 19 high resolution QSO spectra, such that: 
(F) bs = e' T "« = [0.878±0.019,0.696±0.025,0.447±0.031]. 
The uncertainties on these values were estimat ed by binning 
the la errors on r e g presented in Table 5 of ISchave et al] 
(2003) into redshift bins of width Az = 0.5. These values 
have been corrected for absorption from metal lines and 
Lya systems with damping wings. We rescale the optical 
depth in each pixel by a constant factor A to reproduce 
(F 1 )obs at the appropriate redshift. Hence, for a set of spec- 
tra with a total of N pixels: 

JV 

(F) obs = l^e-^=e-^ (3) 
3=1 

where 7] is optical depth in each pixel of the simulated spec- 
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tra. The value of F_i2 required to reproduce {F) b a from 
the simulated spectra is then calculated. 

2.4 The effect of box size and mass resolution on 
the ionization rate 

We firstly vary the box size and resolution of simulations 
of our fiducial model to assess their effect on the value of 
T_i2 . The simulation box must be large enough to include 
long wavelength density perturbations and represent a fair 
sample of the Universe while retaining adequate resolution. 
The left panel of fig.^shows r_i2 for different box sizes at 
fixed mass resolution. As the box size increases, the simu- 
lated optical depths decrease; low density gas responsible for 
Lya absorption falls into the haloes formed in the simula- 
tion, producing larger and emptier voids. A smaller F_i2 is 
thus required to reproduce the observed optical depth. The 
data have almost converged for a box size of 30/i _1 comov- 
ing Mpc, with a difference of less than 2 per cent between 
the T-12 inferred from the 60 - 400 and 30 - 200 models. 

The right panel illustrates the effect of changing the 
mass resolution on T_i2 for fixed box sizes of 15/i _1 and 
30/i~ comoving Mpc. Increasing the mass resolution de- 
creases the optical depth as higher resolution simulations 
resolve smaller haloes. More gas is taken from the low den- 
sity IGM and a smaller value of F_i2 is required to match 
the observed mean flux. Simulations with less than 400 3 
gas particles within a 30/i _1 comoving Mpc box do not 
achieve adequate numerical convergence. The difference be- 
tween F_i2 calculated from the 200 3 and 400 3 simulations 
is on average 7 per cent. Comparison of the simulation with 
200 3 and 400 3 particles within a 15/i _1 comoving Mpc box 
shows that a further increase in resolution reduces F_i2 by 
an extra 5 per cent. 



Unfortunately, 400 3 simulations are too computation- 
ally expensive for the purpose of exploring a large parame- 
ter space. We compromise by taking our fiducial model to 
have 200 3 gas particles within a 15/i _1 comoving Mpc box; 
the value of T_i2 we infer will be systematically high by 
around 8 per cent due to the combined error from box size 
and resolution. Hereafter, all results based on simulations 
with 200 3 particles within a \5li~ 1 comoving Mpc box will 
be reduced by this factor to correct for box size and resolu- 
tion. Note that this correction is likely to be a lower limit 
to the true correction required; our highest resolution sim- 
ulations may not yet have reached numerical convergence. 
We will account for this uncertainty in our error analysis in 
section 15. f I 



3 SCALING RELATIONS FOR T 12 FROM 
SIMULATED ABSORPTION SPECTRA 

It is implicit in the FGPA scaling relation that the den- 
sity and velocity distribution, along with the effective equa- 
tion of state of the Lya absorbers remain unchanged for 
different values of fit,, h, fi m , T and F_i2. Hence, when 
inferring r_i2 from hydrodynamical simulations and subse- 
quently rescaling with different parameters, it is generally 
assumed that: 

r_i2 oc figft 3 T-°- 7 fi-°- 5 (4) 

(e.g. Rauch et al. 1997, McDonald & Miralda-Escude 2001). 
However, it is not clear if this relation is an accurate approx- 
imation. Altering fi m is likely to alter the gas density dis- 
tribution, and thermal broadening may produce a deviation 
from the T -0 ' 7 scaling based on the temperature dependence 
of the recombination coefficient. In addition, the dependence 
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of r_i2 on the polytropic index 7 of the effective equation 
of state T = TbA 7-1 is unclear. Although hydrodynamical 
simulations predict 7 ~ 1.6 at low density, i f He 11 reion- 
ization occurs at z ~ 3 l|Reimers et al.lll997l : ISchave et alJ 
l200d: iBernardi et all 12003) " the effective equation of state 
for the low density gas may become isothermal. Radiative 
transfer effects may also raise the temperature of the gas 
and blur the monotonic relation between temperature and 
densi ty jAbel fc Haehneltl Il999l iBolton. Meiksin fe White! 
12004) . Finally, changing ag will affect T_i2 by changing the 
density distribution of the gas. 

We test the validity of equation 01 using our suite of 
hydrodynamical simulations. The motivation for this is two- 
fold: first, it is important to understand the physical pro- 
cesses which alter the value of T_i2 we infer from simula- 
tions; second, we wish to determine how F_i2 scales with 
the input parameters, allowing us to place joint constraints 
upon them. We do not discuss the scaling with f2b and h, 
since our analysis revealed that the FGPA scaling relation 
holds extremely well for these parameters. 

3.1 The dependence of T 12 on temperature 

The left panel of fig. [5] shows our test of the FGPA scal- 
ing relation for the gas temperature at mean density using 
five simulations with photoheating rates [0.3,1.0,2.0,3.3,4.0] 
times the HM96 values. We expect the ionization rate to de- 
crease as the HM96 photoheating rate is raised; the hydro- 
gen recombination rate scales as T~ 0,7 , so for larger tem- 
peratures there is less neutral hydrogen. We find a mod- 
est deviation from the T -0 ' 7 scaling (shown as the dotted 
curves) at low temperatures, probably due to thermal broad- 
ening of absorption features, an effect which is strongest 
at z = 4. The fiducial model with its higher gas tempera- 
tures will have relatively broader absorption features, pro- 
ducing more Lya absorption and thus requiring a larger 
value of F_i2 than if thermal broadening were absent. The 
FGPA scaling relation for the fiducial model will slightly 
over-predict the value of F_i2 at lower temperatures. The 
same process also occurs for gas temperatures higher than 
the fiducial model; the absorbers are thermally broadened 
to a greater extent so the scaling relation under-predicts 
the ionization rate. This effect is much more apparent at 
higher redshifts when the spatial density of Lya absorbers 
is greater. Note that the temperature dependence of the gas 
pressure may also contribute to the deviation from the T~ ' 7 
scaling. We perform a least squares fit on the data to obtain 
the following scaling for T_i2 with temperature: 

F_ 12 ocT xl(z) (5) 

where x 1 (z) = [-0.68, -0.62, -0.59] at z = [2, 3, 4]. We shall 
adopt these scaling relations (shown as solid curves) from 
here onwards. 

3.2 The dependence of T_i2 on 7 

The dependence of the hydrogen ionization rate on the poly- 
tropic index 7 is shown in the right panel of fig.|5| The value 
of 7 is varied by artificially rescaling the effective equation 
of state of our existing fiducial simulation data by pivot- 
ing the temperature-density relation around the mean gas 



density JViel. Haehnelt fc SpringejEoolh . Note that the ac- 
tual slope of the effective equation of state predicted by the 
numerical simulations is 7 ~ 1.6. This will not account for 
any dynamical effects the different temperatures have on 
the gas distribution, but it will model thermal broadening 
and the change in neutral hydrogen density correctly, thus 
providing a reasonable approximation. Note that the data 
point with the largest 7 is for our simulation of the fiducial 
model with no scaling of 7. It is clear from fig. that the 
inferred ionization rate has a non-negligible dependence on 
7. At z = 3 and z — 4 T_i2 increases with 7. A flatter 
effective equation of state produces hotter gas at densities 
less than the cosmic mean with a lower neutral hydrogen 
fraction. As most of the Lya absorption is produced by this 
gas, decreasing 7 will decrease the amount of absorption and 
hence lower the inferred ionization rate. The dependence of 
r_i2 on 7 is much stronger at z — 4 compared to z = 3. 
However, at z — 2, the trend is opposite; the ionization rate 
decreases with increasing 7, albeit with a weak dependence. 
This can be explained by considering the density evolution 
of the Lya absorbers. With decreasing redshift, the mean 
gas density decreases. By z = 2, the optical depth in re- 
gions of average density and below has become so low that 
these contribute very little to the overall effective optical 
depth. The absorption is almost entirely due to overdense 
regions. Decreasing the polytropic index will lower the tem- 
perature of this gas, increasing the neutral hydrogen fraction 
and slightly increasing the effective optical depth. The in- 
ferred ionization rate increases slightly with decreasing 7 as 
a consequence. In a similar manner to before, we perform a 
least squares fit on the data to obtain a scaling relation for 
the ionization rate with the polytropic index 7: 

r_ 12 oc 7^ (6) 

where x 2 (z) = [-0.04,0.34,0.55] at z = [2,3,4]. These scal- 
ing relations are shown in fig. [5] as the solid curves. 

3.3 The dependence of T 12 on Q m 

The FGPA scaling relation assumes that the ionization rate 
per atom scales as Q^ ' 5 . For this to be true, the density 
distribution of the neutral hydrogen gas must remain the 
same for different values of Q m . We test this scaling relation 
using five simulations with Q m — [0.17, 0.26, 0.4, 0.7, 1.0] and 
show the results in the left panel of fig. El As before, the plot- 
ted values are rescaled to match the fiducial temperatures, 
7 = 1.3, corrected for box size and resolution and are shown 
with the best fitting curves to the data. The f2m°' 5 scaling 
around the fiducial model is indicated by the dotted lines. 
This becomes a poorer fit to the data as f2 m is increased, 
and overestimates the value of F_i2 substantially. It is worth 
stressing that all the models shown in the left panel of fig. [3] 
have the same normalisation at a scale of 8/i _1 Mpc for the 
dark matter linear power spectrum (as = 0.85). In order to 
understand why the discrepancy with the expected trend is 
so big we performed a new run for the Einstein-de Sitter 
(M4) model and normalised the power spectrum to have the 
same fluctuation amplitude as the fiducial model at the scale 
of 30 kms" 1 , roughly corresponding to the Jeans length. The 
filled points in the left panel of fig.|3]indicate the new values 
of r_i2 we infer; one can see that these values are in better 
agreement with the theoretical curve, showing that F_i2 is 
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Figure 2. Left: The dependence of T_i2 on the gas temperature at mean density To at three different redshifts as indicated on the plot. 
Results are corrected for box size and resolution as discussed in section 12.41 Solid curves show a least squares fit and the dotted curves shows 
the T~ 0,7 scaling due to the temperature dependence of the recombination coefficient. The filled symbols show the fiducial temperatures 
and their uncertainties. Right: The dependence of T_i2 on the index of the temperature density relation 7. Results are corrected for box 
size and resolution. The filled circle shows our fiducial value of 7 and its uncertainty. 




Figure 3. Left: The dependence of the estimated T_i2 on O m for three different redshifts as indicated on the plot. The solid curves are 
a least-square fit. The dotted lines show the O^ ' 5 scaling. Filled points are obtained for a model with f2 m = 1.0 and the the same r.m.s. 
fluctuation amplitude at a scale of 30 kms -1 as our fiducial model with Q m = 0.26. As before the values are corrected for box size and 
resolution. The filled circle shows our fiducial value of f2 m = 0.26 and its uncertainty. Right: The dependence of T_i2 on ag. The dotted 
curve shows the result of Tytler et al. at z = 1.9, assuming {F) ^ s = 0.882 ± 0.01. Results are corrected for box size and resolution. The 
filled circle shows our fiducial value of erg and its uncertainty. 
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Figure 4. The PDF of the volume weighted gas density as a 
fraction of the mean at z = 2 for the fiducial Q m = 0.26 (solid 
line) and M4 f2 m = 1.0 (dotted line) models. The fiducial model 
has a gas density distribution which peaks at a greater amplitude 
closer to the mean density, producing more Lya absorption at 
fixed r_i2. 
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Figure 5. The dependence of T_i2 on the effective optical depth 
at different redshifts for our fiducial model with our fiducial values 
of To and 7. The symbols indicate the values of r ef f quoted by 
various authors. The T g from Tytler et al. (2004) is for z = 1.9. 
The filled circles indicate our fiducial parameter values, with the 
estimated uncertainty shown by the horizontal line. 



more sensitive to the r.m.s. fluctuations on a fixed velocity 
scale. This suggests that the r.m.s. fluctuation amplitude at 
the Jeans scale is more relevant. 

Fig. [I] shows the probability density function (PDF) of 
the volume weighted gas density as a fraction of the mean 
at z — 2 for our fiducial model and the Einstein-de Sit- 
ter (M4) model. The density distribution of the Einstein-de 
Sitter model peaks at lower density, explaining the much 
lower P_i2 required to match the observed mean flux. In the 
Einstein-de Sitter model a greater proportion of the gas has 
fallen into haloes. Different gas density distributions are not 
accounted for in the FGPA scaling relation. Comparisons 
of parameters such as T_i2 and Qbh 2 inferred from simu- 
lations with different matter distributions using the FGPA 
gives misleading results. A least squares fit to our data yields 

r_ 12 cx^ (z \ (7) 

where xs(z) = [-1.00, -1.04, -1.16] at z = [2, 3, 4]. 

3.4 The dependence of T 12 on as 

The magnitude of as, the mass fluctuation amplitude within 
8/i _1 Mpc spheres, also alters the required F_i2 . For small 
as, there is less small scale power and the collapse of struc- 
tures is less advanced. This produces a narrower volume 
weighted density PDF with a larger average density. This 
decreases the mean flux of the simulated spectra if the ion- 



ization rate is kept fixed. A larger F_i2 is thus required to 
reproduce the observed effective optical depth. 

The right panel of fig. [3] shows F_i 2 for as = 
[0.5,0 7,0.85,1.0,1.21 at redshifts z = [2. 3. 41. FlVtler et alJ 
(2004) present a similar plot in their fig. 16, where they 
construct curves relating as and 7912 = r_i2/1.329, nor- 
mali sed to the predicted ioni zation rate at z = 1.9 from 
iMadau. Haardt fc Reesl lll999l) . The central curve in their 
plot gives the T_i2 and as they require to reproduce 
(F} hs = 0.882 ± 0.01. We find a weaker dependence of 
F_i2 on (78. We use similar values for f2 m , h, fib, the UV 
heating rate, and the slightly lower value of r c ff , which Tytler 
ct al. use, will change our results only by a small amount 
(see section l3~5t . The small difference in redshift should also 
not alter the results significantly, and the resolution and 
box size of our simulations will also not account for the 
difference; increasing these systematically lowers F_i2. The 
most likely reason for this discrepancy are differences in the 
modelli ng techniques for ou r respective cosmological simu- 
lations (lO'Shea et al]l2003D. Tytler e t al. use the Eulerian 
code ENZO dNorman fc Brvanlll999l) . This models the gas 
component using a piecewise parabolic method, whereas we 
use a Lagrangian SPH code. A least squares fit to our results 
gives: 

F_ 12 oc a^ (z) (8) 
where x A (z) = [-0.90,-0.99,-1.26] at z = [2,3,4]. 
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3.5 The dependence of T_i2 on r e g 

To determine the intensity and evolution of F_i2 using the 
FGPA we must assume a value for the effective optical 
depth, T c ff = — ln(F) b s . Measurements of the effective op- 
tical depths have been made by a number of authors (e.g 
Press, Rybicki & Schneider 1992, Kim et al. 2001, Songalia 
& Cowie 2002, Bernardi et al. 2003, Schaye et al. 2003). 
The values obtained from low-resolution, low-S/N data by 
Press, Rybicki & Schneider and Bernardi et al. appear to 
be systematically higher than those obtained by the other 
groups, mainly from high-resolution data (see Viel, Haehnelt 
& Springel 2004 for a discussion). The differences are pri- 
marily due to the systematic uncertainties in the continuum 
fits. Tytler et al. make an analysis of continuum fitting errors 
by comparing the true and fitted continua for sets of artifi- 
cial spectra. They find that they can obtain a fit which is on 
average good to within 1-2 per cent, except for very low S/N 
spectra. Some agreement is emerging that the lower values 
from the higher quality data are more appropriate. Statis- 
tical errors from the presence of metal absorption features 
and high column density absorber s also contribute to the 
uncertainty in r e ff llViel et al.ll2u0i) . 

Fig. shows the value of F_i2 we infer from our fidu- 
cial model at z = [2, 3, 4] by varying the effective optical 
depth, and hence mean observed flux, to which the simu- 
lated optical depths are scaled to match. At z — 2, a change 
in r e g of 0.05, which is consistent with the variation in es- 
timates, results in almost a factor of two difference in the 
required hydrogen ionization rate. An accurate estimate of 
T e ff and the errors associated with the measurement are re- 
quired to determine F_i2; at z = 2 even a small change 
in the assumed r e s will change the inferred ionization rate 
appreciably. The best fitting curve to the data is: 

r-12 oc r3« (9) 
where x 5 (z) = [-1.44, -1.61, -1.68] at z = [2, 3, 4]. 

4 COMPARISON WITH PREVIOUS 
ESTIMATES FROM SIMULATIONS 

4.1 Hydrodynamical codes 

In this section we shall make a direct comparison between 
our data and the results of other authors. Firstly, we con- 
sider two simulations which are similar to the ACDM and 
SCDM models analysed by Rauch et al. (1997). The simu- 
lations have identical cosmological parameters, box size and 
resolution. Our ACDM model differs from that analysed by 
Rauch et al., since we use an SPH code rather than an Eule- 
rian grid code. We try to approximate the different temper- 
atures of the models by increasing the photoheating rate to 
3.3 times the HM96 value for our ACDM simulation. This 
results in similar temperatures as those given by Rauch et 
al. at z = 2 (this is only redshift for which the temperature 
was published). The optical depths of our artificial spectra 
are scaled to match the mean flux used by Rauch et al., cor- 
rected for continuum fitting: (F) obs = [0.848, 0.670, 0.414] 
and [0.846, 0.655, 0.383] at z = [2, 3, 4] for the ACDM and 
SCDM models respectively. Mimicking Rauch et al., we 
rescale F_i2 by a factor (fl^h 2 ) 2 , setting Q,^h 2 — 0.024. Un- 
like in the rest of the paper, we have not corrected our data 
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Figure 6. Comparison of the ionization rate inferred from 
GADGET-2 to that inferred from the simulations used in Rauch 
et al. (1997). Diamonds and triangles represent the ACDM and 
SCDM models, respectively. The GADGET-2 results are shown 
with solid lines while the results from Rauch et al. are shown with 
dotted lines. 

for numerical effects or rescaled the temperatures to match 
our fiducial values of To and 7. 

Fig. HJ shows the redshift dependence of T_i2 inferred 
from the ACDM and SCDM models as diamonds and tri- 
angles respectively. The solid lines give our data and the 
Rauch et al. data are plotted as dotted lines. Rauch et al. 
conclude the large difference between F_i2 obtained using 
the ACDM and SCDM models is primarily due to differ- 
ent gas temperatures in the respective simulations. The gas 
temperature at mean density in their SCDM model is about 
a half of that in the ACDM model at z = 2. We suggest the 
values of Q m , erg and 7 used in each of the simulations will 
also alter F_i2. The different amounts by which the optical 
depths of the artificial spectra are scaled will also produce 
a small difference in F_i2, but this change will be negligible 
compared to the effect of the other parameters. In our sim- 
ulations the smaller temperature of the SCDM simulation 
counteracts the effect of different fi m and as and we find 
rather similar values for both simulations. The values for 
our ACDM simulation are larger than those in Rauch et al. 
while those for the SCDM simulation are smaller. Unfortu- 
nately the only information about the thermal state of the 
gas provided by Rauch et al. is at 2 = 2 which makes a ac- 
curate comparison difficult. Although we attempt to match 
the gas temperature in our models at 2 = 2 to those given by 
Rauch et al., this agreement is probably not good at higher 
redshifts, particularly for the ACDM model. The gas tem- 
peratures of the Rauch et al. models result from different 
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reionization histories instead of artificially increased photo- 
heating rates. The UV background of their ACDM model 
ionizes He n at z ~ 3, which may lower the polytropic in- 
dex 7. The 7 in the Rauch et al. ACDM simulation may thus 
have been lower than that in our simulations. This would 
make the discrepancy smaller. 

ISchave et alJ J2003) utilise a hydrodynamical simula- 
tion based on a modified version of t he SPH code HYDRA 
iCouchman. Thomas fc Pearcdll995ll to examine the metal- 
licity of the IGM. Their model consists of a periodic box of 
12/i -1 comoving Mpc with 256 3 dark matter particles and 
256 3 gas particles. The parameters they adopt are fi m = 0.3, 
tt h h 2 = 0.019, h = 0.65, a$ = 0.9 and n = 1.0 with a 
UV background designed t o reproduce the effe ctive equa- 
tion of state constrained bv lSchave et alJ i2000h . The opti- 
cal depths of the artificial Lya spectra are scaled to match 
the same mean flux as ours. They find the P-i2 required 
to match the mean flux of the Lya forest at z = [2, 3, 4] 
are F_i2 = [0.87,0.54,0.36]. Rescaling our fiducial model 
values of r_i2 = [1.29,0.86,0.97] to match the cosmologi- 
cal parameters used by Schaye et al. reduces our values by 
around 40 per cent to T_i 2 = [0.80,0.53,0.57]. We remove 
the correction for box size and resolution from these val- 
ues to allow a fair comparison. These are consistent with 
ISchave et alJ (Hp03), allowing for the unknown differences 
between the exact thermal state of the low density gas in 
the two simulations. Theuns, Schaye & Haehnelt (2000) also 
run simulations using HYDRA to study the broadening of 
Lya absorption lines in QSO spectra. Their reference model 
consists of 64 3 dark matter particles and 64 3 gas particles 
within a periodic box of 2.5h~ 1 comoving Mpc, with the 
same cosmological parameters used by Schaye et al. and 
the UV background specified by Haardt & Madau (1996). 
This produces a temperature of To ~ 11750 K at mean 
cosmic density. To reproduce an effective optical depth of 
i~cff = 0.33 at z = 3 they require P_i2 = 0.87, based on the 
value of A quoted in their paper. Scaling our fiducial data 
to match their parameters reduces our value by 10 per cent 
to r_i2 = 0.79. Although the mass resolution of our simu- 
lations is comparable, we obtain a smaller value because of 
our larger simulation box size. 

Most recentlv. lTvtler et alJ l)2004h made a detailed anal- 
ysis of the Lya forest at z = 1.9, using a set of 77 moderate 
resolution QSO spectra and compared it to a large hydro- 
dynamical simulation run with the Eulerian code ENZO. 
Including uncertainties on the value of erg, fib and r e s , they 
find a joint constraint of T_i2 = 1.44 ± 0.36. They use a 
hydrodynamical simulation on a 1024 3 grid in a 54.528/I." 1 
comoving Mpc box, with fi m =0.27, Q b h 2 = 0.022, h = 0.71, 
as = 0-9 and n = 1.0 a n d a U V background specified 
by iMadau. Haardt fc Reesl (^999) with an artificially in- 
creased He 11 photo-heating rate. Our fiducial result at 
z — 2 rescaled to match the cosmological parameters and 
i~cff = 0.126 at z = 1.9 assumed by Tytler et al. gives 
r_i2 = 1.17. This value is in reasonable agreement with 
the Tytler et al. result, again allowing for differences be- 
tween the exact gas temperatures and numerical method. 
Note, however that Tytler et al. found a different scaling of 
r_i2 with as (see section l3~4l for details). 



Table 3. Fiducial parameter values and estimates of their uncer- 
tainty. The values listed for gas temperature at mean density To 
and the effective optical depth r e g are for z = [2, 3, 4]. 



Parameter Fiducial values and uncertainties 



T [11200, 17800, 12500]±5000 K 

£V 0.26±0.04 

r off [0.130±0.021, 0.362±0.036, 0.805±0.070] 

7 1.3±0.3 

U h h 2 0.024±0.001 

a 8 0.85±0.05 

h 0.72±0.04 



4.2 Pseudo-hydrodynamical codes 

Several authors have also considered this problem using 
pseudo-hydrodynamical codes (e.g Gnedin & Hui, Gnedin 
1998, Meiksin & White 2003, 2004). These approximate 
methods sacrifice accuracy for speed, allowing a much larger 
parameter space to be explored. Providing the accuracy of 
these approximations and the parameter ranges they model 
effectively are well controlled, they c an provide an alter- 
native method for modelling the IGM. iGnedinl jl998T) con- 
cludes that observational constraints on the column density 
distribution of the Lya forest are consistent with the inten- 
sity of the UV b ackground remaining cons tant in the redshift 
range z ~ 2 - 4. IMeiksin fc White! i2004) . in a study of UV 
background fluctuations, find the mean flux of the IGM is 
reproduced with r_i 2 = [0.88^;^, 0.76± ;^] at z = [3,4]. 
They adopt fi m =0.3, fi b /i 2 = 0.020, h = 0.70, a 8 = 1.0 
and n = 1.0, and impose and e ffective equation o f state 
with T = 20000 K and 7 = 1.5. IMeiksin fc White! (1200 ill 
find that the spectral properties of the Lya forest obtained 
from their code are reproduced to within 10 per cent of those 
found using fully hydrodynamical simulations. These results 
can thus be considered to be consistent with ours. 



5 THE METAGALACTIC HYDROGEN 
IONIZATION RATE 

5.1 The error budget 

We now proceed to place joint constraints on the metagalac- 
tic hydrogen ionization rate using independent estimates 
for several cosmological and astrophysical parameters. We 
adopt the values listed in Table |3] The values of To and 7 
are based on the results of Schaye et al. (2000) and we as- 
sume the Lya effective optical depth advocated by Schaye et 
al. (2003). The cosmological parameters are consistent with 
the results of Spergel et al. (2003), based on their analy- 
sis of the cosmic microwave background anisotropy. Using 
our scaling relations, we estimate the error on the values 
of r_i2 from our fiducial simulation, corrected for box size 
and resolution. The total error is calculated by adding the 
individual contributions in quadrature, with the percentage 
error on T_i2 contributed by each parameter listed in Ta- 
ble^] We make the conservative assumption of an additional 
numerical uncertainty of 10 per cent to take into account 
numerical convergence and the modest differences with nu- 
merical simulations of other authors. The total percentage 
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Figure 7. Comparison of our best estimate for T_i2 at z = [2,3,4] with constraints from observations. Our data are plotted with filled 
diamonds, and the grey shaded area shows the error bounds. The filled squares and triangles show the estim ated contribution to the 
metaealactic ionization rate from QSOs+galaxies and QSOs alone, based on estimates from the updated model of lMadau. Haardt &; Reesl 

1999) including UV photons from re-processing by the IGM. The hatched box gives the constraint on T_i2 from the proximity effect 
i Scot^e^alJl2000|) and the estimate from Lyman-break galaxies assuming a global spectral index of a = 1.8 is plotted with a filled circle 

Stcid el et al .1200 J) . ' me data po int has been offset from z = 3 for clarity, and the upper error limit is not shown. The star shows the best 
estimate of T_i2 at z = 1.9 from lTvtler et all l2004h . including their errors from the uncertainty in T e g, erg and H^. 



Table 4. Percentage error budget for r_i2 from estimates of 
various cosmological and astrophysical parameters, listed approx- 
imately in order of importance. The total error is obtained by 
adding the individual errors in quadrature. 



Parameter 


2 = 2.00 


z = 3.00 


z = 4.00 


T 


+ 50 of 
-22 /o 


+23 Of 
-14 /o 


+35 of 
-18 /o 




-13 /o 


+19<w 
-14 /c 


+21 of 
-15 /o 


T e K 


+29 of 


+18 of 
-14 /o 


-13 /o 


Numerical 


±10% 


±10% 


±10% 


7 


±1% 


+7<w 
-9 /o 


tl 2 3% 


n h h 2 


+9of 
-8 /o 








+Gof 
-5 % 


±6% 




h 


±6% 


±6% 


±6% 


Total 


-36 /o 


+39 of 
-30 /o 


+49 of 
-34 /o 



errors on T_i2 are [±§§,±^,±34] per cent at z — [2,3,4], 
with the largest contribution coming from the uncertainty 
in the gas temperature. There is also a large contribution 
from the uncertainty in the effective optical depth, espe- 



cially at z = 2. Interestingly, the uncertainty in fi m also 
gives a substantial contribution to the total error budget. 
This is primarily due to the sensitivity of F_i2 on f2 m due 
to changes in the r.m.s. fluctuations of the gas density at the 
Jeans scale. Other uncertainties are less important, in par- 
ticular Q.hh 2 , as and h. This contrasts with the conclusion of 
Tytler et al., who find as is the biggest contribution to the 
uncertainty on r_i2. However they find a stronger depen- 
dence of F_i2 on as than we do, and they do not investigate 
the contribution of uncertainties in Q m or T to this error. 
Note also that inhomogeneities in the UV background may 
increase estimates of F_i2 by up to twenty per cent if the 
UV background is dominated by emission of QSOs (Gnedin 
& Hamilton 2002; Croft 2004; Meiksin & White 2004). We 
will conclude in section 5.2. that QSOs contribute less than 
half of the ionizing UV background. We will therefore ne- 
glect a possible error due to inhomogeneities in the UV back- 
ground. Using our fiducial model, we find the Lya effective 
optical depth of the IGM at z = [2, 3, 4] is reproduced by 

r_i2 = [i.29±8:t§,o.86±8:iS,o.97±8:g]. 



5.2 Comparison with other observational 
estimates 

We compare our results with other observational estimates 
of r_i2 in fig. [7| Our estimates are plotted as filled dia- 
monds, with the error bound shown by the dark grey shaded 
region. The filled triangles show the metagalactic hydro- 
gen ionization rate computed from the updated UV back- 
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ground model of iMadau. Haardt fe Reesl <|l999!) , hereafter 
MHR99 (Francesco Haardt, private communication). The 
values quoted here are smaller by more than a factor of 
two than those presented in the work of HM96, primarily 
due to different assumptions about the emissivity of the 
ionizing sources, and to a lesser extent upon the assumed 
cosmology. The rates are based on the QSO optical lumi- 
nosity function of Boyle et al. (2000), with a redshift evo- 
lution similar to that of MHR99, aside from a steeper evo- 
lution for z > 3 to fit data from the SDSS survey (Fan 
et al. 2001). These values are computed assuming the cur- 
rent cosmological concordance model, and also include the 
contribution of diffuse radiation from re-processing by the 
IGM (e.g. HM96). Our data appear to be inconsistent with 
the IGM being kept ionized by QSOs at z < 4. For com- 
parison, the filled squares give F_i2 calculated using the 
QSO rate above, plus an additional source in the form of 
Young Star Forming Galaxies (YSFGs). The YSFG rates 
are based on the models of Bruzual & Chariot (1993), as- 
suming a Salpeter IMF, a constant star formation rate, 
an age of 500 Myr, Z = 0.2Z@ and a star formation his- 
tory as in MHR99. The very uncertain escape fraction is 
assumed to be / csc = 0.1. Note tha t at low redshift the 
observed value is less than half this (fLeitherer et aljfl995t 
iFernandez-Soto. Lanzetta fc Chenll2003F) . 

These ionization rates are in good agreement with 
our results, suggesting that a substantial contribution from 
galaxies appears to be required at all redshifts. Note, how- 
ever that there are still substantial uncertainties in the es- 
timates of the ionization rate due to QSOs and galaxies. 
Hunt et al. (2004), Croom et al. (2004) and Barger et al. 
(2004) all suggest that the faint end slope of the luminosity 
function of optically and UV bright QSOs is shallower than 
the extrapolation from the Boyle et al. data, which may fur- 
ther lower the expected contribution by QSOs. On the other 
hand Scott et al. (2004) find that the UV spectrum of low 
luminosity QSOs is substantially harder than that of bright 
QSOs which may somewhat increase their contribution to 
the UV background. The main uncertainty in the estimate 
for the ionization rate due to the UV emission of galax- 
ies is the escape fraction of ionizing photons. However, the 
improved estimates for the QSO emissivity are so low that 
there should now be little doubt that the UV background 
cannot be dominated by QSOs over the full redshift range 
z < 2 < 4. Aguirre et al. (2004) come to similar conclusions 
in their study of the Si iv /C iv ratio of the associated metal 
absorption in Lya forest spectra. 

Constrai nts on r_i2 from the proximity effect mea- 
surements of lScott et al] (120001) are shown by the hatched 
box, based on a set of 99 moderate resolution QSO spectra. 
Similar results have been obtained by Cooke et al. (1997). 
Allowance has been made for variation in the spectral in- 
dices of QSOs and differences in systemic redshifts measured 
from the [O in] or Mg n line and Lya emission. Our results 
are consistent with the lower end of these estimates. The 
filled circle gives the F_i2 we calculate assuming a spec- 
tral index of a = 1.8 using the metagalactic ionizing radia- 
tion intensit y, J^(912 A), inferr ed from Lyman-break galax- 
ies (LBGs) Isteidel et alll200ll) . Steidel et al. quote J„(912 
A) ~ 1.2 ±0.3 x 10" 21 ergs" 1 cm" 2 Hz" 1 sr" 1 at z ~ 3 based 
on observations of the bluest 25 per cent of LBGs. We offset 
the data point from z = 3 for clarity, and the upper error 



limit extends beyond the range of the plot. We estimate the 
maximum error on this result to be a factor of 4 to allow for 
the possibility that the 75 per cent of LBGs which are redder 
may contribute negligibly to the UV background. Note that 
Giallongo et al. (2002) found an upper limit to the escape 
fraction from two bright LBGs which is about a factor four 
lower than that measured by Steidel et al.. There are further 
uncertainties due to the mean free path of ionizing photons 
assumed by Steidel et al. and the spectral index that we 
assumed. 

Finally we als o plot the value of F_i2 inferred by 
iTvtler et al ] J2004I) in their study of the Lya forest at 
z = 1.9. The error bars include the uncertainties in T e g, 
us and Qb , based on scaling from their hydrodynamical sim- 
ulations. 



6 CONCLUSIONS 

We have examined the dependence of the ionization rate per 
hydrogen atom on cosmological and astrophysical parame- 
ters using hydrodynamical simulations. This is the first time 
all these parameters have been studied using fully hydrody- 
namical SPH simulations with high resolution in large boxes. 
We find that the gas temperature, the effective optical depth 
Tefr and the mass fraction Sl m have a marked effect upon the 
hydrogen ionization rate, in addition to the baryon fraction 
Qh, ""8, 7 and the Hubble constant. We find that F_i2 scales 
around our fiducial model as 

F_i2 cx Qfch 3 T clW i* 2W SQWaZ* (z) T3i w , (10) 

where we tabulate the indices of equation (10) in Table |S] 
We stress that this scaling relation is likely to be somewhat 
model dependent and should not be applied to models with 
parameters very different from our fiducial model without 
further checks. However, it does provide a clear picture of the 
degeneracies which exist between F_i2 and other parameters 
within our simulations, and is independent of assumptions 
about the gas distribution, effective equation of state and 
ionized gas fraction of the low density IGM. In particular, 
the Lya optical depth at fixed r.m.s. fluctuation amplitude 
as is more strongly dependent on f2 m than the FGPA sug- 
gests due to changes in the gas density distribution. Thermal 
broadening also produces a deviation from the FGPA scaling 
for temperature, although this change is much less dramatic. 
Consequently, we urge caution when using the FGPA scal- 
ing relation, especially if comparing models with differing 
ttm. We find good agreement with estimates of F_i2 from 
the Lya forest opacity by other authors if differences in the 
assumed cosmological, astrophysical and numerical param- 
eters are taken into account. 

In recent years, compelling observational evidence has 
led to the acceptance of a standard cosmological model 
which is flat, has low matter density and a substantial con- 
tribution of dark energy to the total energy density. Within 
this 'concordance' cosmological model the current genera- 
tion of hydrodynamical simulations predict values for the 
metagalactic hydrogen ionization rate, required to reproduce 
the effective Lya optical depth of the IGM in the range 
z = 2 — 4, which are about a factor of four larger than 
those in an Einstein-de Sitter model with the same r.m.s. 
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Table 5. The redshift dependent indices from our scaling relation 
of r_i2 with several cosmological and astrophysical parameters. 





T 


7 


a2m 


CT8 




2 


xi(z) 


X2{z) 


x 3 (z) 


xa(z) 


X5(z) 


2 


-0.68 


-0.04 


-1.00 


-0.90 


-1.44 


3 


-0.62 


0.34 


-1.04 


-0.99 


-1.61 


4 


-0.59 


0.55 


-1.16 


-1.26 


-1.68 



density fluctuation amplitude erg. The ionization rates esti- 
mated from the Lya forest opacity are more than a factor 
two larger than estimates from the integrated flux of op- 
tically/UV bright observed QSOs alone. This discrepancy 
increases with increasing redshift. We confirm the findings 
of Tytler et al. (2004) at z ~ 1.9 that the estimated ioniza- 
tion rates from simulations of a ACDM concordance model 
are in reasonable agreement with the estimates of the inte- 
grated ionizing flux from observed QSOs plus a significant 
contribution from galaxies as in the model of MHR99. A 
model where the UV background is dominated by emission 
from QSOs appears to be inconsistent with the Lya forest 
data in the complete redshift range 2 < z < 4. The estimates 
of the ionization rate are also in agreement with the lower 
end of the range of estimates from the proximity effect. 

There are still many issues which need to be resolved. 
In particular, better constraints on the thermal history of 
the IGM are required. Additional physics such as radiative 
transfer, galactic feedback and metal enrichment may need 
to be incorporated into simulations in a more realistic fash- 
ion. However, our estimate of the errors on F_i2 most likely 
accounts for the modest changes expected from these pro- 
cesses, suggesting we have obtained a consistent constraint 
on the metagalactic hydrogen ionization rate. 
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